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Abstract 

In the present work, we discuss a procedure for constructing sparse binary and ternary 
matrices from existing two binary sensing matrices. The matrices that we construct have 
several attractive properties such as smaller density, which supports algorithms with low 
computational complexity. As an application of our method, we show that a CS matrix of 
general row size different from p,p 2 ,pq (for different primes p,q) can be constructed. 
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1. Introduction 

Recent developments at the intersection of algebra and optimization theory, by the name 
of Compressed Sensing (CS), aim at providing sparse descriptions to linear systems. These 
developments are found to have tremendous potential for several applications [7, 13, 21, 26, 
27]. Sparse representations of a function are a powerful analytic tool in many application 
areas such as image/signal processing and numerical computation [8], to name a few. The 
need for the sparse representation arises from the fact that several real life applications 
demand the representation of data in terms of as few basis (frame) type elements as possible. 
The elements or the columns of the associated matrix <I> are called atoms and the matrix so 
generated by them is called the dictionary. The developments of CS Theory depend typically 
on sparsity and incoherence [8] [25]. Sparsity expresses the idea that the information rate 
of a continuous time data may be much smaller than suggested by its bandwidth or that a 
discrete-time data depends on a number of degrees of freedom which is comparably much 
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smaller than its (finite) length. On the other hand, incoherence extends the duality between 
the time and frequency contents of data. 

Since the sparsity of the CS matrix is key to minimizing the computational complexity 
associated with the matrix-vector multiplication, it is desirable that a CS matrix has smaller 
density. Here, by density, one refers to the ratio of number of nonzero entries to the total 
number of entries of the matrix. The sparse CS matrix may contribute to fast processing 
with low computational complexity in Compressed Sensing [21]. 

In the recent literature on CS [1, 6, 16, 19, 22], deterministic construction of CS matrices 
has gained momentum. R. Devore [15] has constructed deterministic binary sensing matrix 
of size p 2 xp r+1 , where p is prime or prime power. The density of this matrix is K S. Li et. al. 
[23] have generalized Devore’s work, constructing binary sensing matrix of size [P\q x q clyG \ 
where q is any prime power and V is the set of all rational points on algebraic curve X over 
finite field F g . The density of this matrix is A A. Arnini et. al. [4] have constructed binary 
sensing matrices using OOC codes. The density of this matrix is A, where m is row size and 
A is the number of ones in each column. In all these constructions row sizes of associated 
matrices are given by some particular family of numbers. 

CS matrices of general size will not impose any restriction on the data to be dealt with, 
when used in applications such as data compression or classification. The main contribution 
of the present work may be summarized as follows: 

• Construction of sparse binary and ternary matrices from existing binary matrices. 

• As an application of our construction methodology, we show that, given m A p,p 2 ,pq 
(for different primes p, q), it is possible to construct sparse binary and ternary matrices 
of row size m. 

Consequently, we believe that this work is an attempt towards constructing CS matrices of 
general size. 

The paper is organized in several sections. In section 2, we provide basics of CS theory. 
While in sections 3 and 4, we discuss in detail the construction of CS matrices of low density 
and an application of proposed method respectively. In the last two sections, we present 
respectively the construction of ternary matrices and concluding remarks. 

2. Basics of Compressed Sensing: 

The objective of compressed sensing is to recover x = (aq, x 2 ..., xm) t G M A/ from a few 
of its linear measurements y G through a stable and efficient reconstruction process via 
the concept of sparsity. From the measurement vector y and the sensing mechanism, one 
gets a system y = 4>a:, where is an m x M (m < M) measurement matrix. 

Denoting <3>£ as the k lLi row of <L, one may rewite the /A 11 component in y as y\~ = 
(Tjh x) , k — 1, 2,... m. Here ($£, x) represents the inner-product between and x. That 
is, the object x to be acquired is correlated with the waveform This is a standard setup in 
several applications [17] [11]. For example, if the sensing waveforms are Dirac delta functions, 
then y is a vector of sampled values of x in time or space domain. If the sensing waveforms 
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are indicator functions of pixels, then y is the image data typically collected by sensors 
in a digital camera. If the sensing waveforms are sinusoids, then y is a vector of Fourier 
coefficients and this modality is used in the magnetic resonance imaging. Nevertheless, 
if the sensing waveforms have 0 and 1 (or 0 and ±1) as elements, then the associated 
matrix (referred conventionally to as a sensing matrix) can have potential application for 
multiplier-less dimensionality reduction. An excellent overview of Compressed Sensing and 
the applicability of several sensing matrices may be seen in [11], 

Given the pair (y,<E>), the problem of recovering x can be formulated as finding the 
sparsest solution (that is, the solution containing most number of zero entries) of linear 
system of equations y = <hx. Sparsity is measured by ||.|| 0 norm and ||:r|| Q = | {j : x 3 ^ 0}|, 
the number of non-zero entries in x. Now finding the sparsest solution can be formulated as 
the following minimization problem (generally denoted as Pq problem): 

Po : nrin ||x||o subject to &x = y. (1) 

X 

This P 0 problem is a combinatorial minimization problem and is known to be NP-hard 
[7]. One may use greedy methods and convex relaxation of P 0 problem to recover A—sparse 
signals (that is, ||cc||o = A: or signals with k number of nonzero entries). The convex relaxation 
of Po problem can be posed as P\ problem [9, 10], which is defined as follows: 

Pi : min||x||i subject to <f>x = y. (2) 

X 

The Orthogonal Matching Pursuit (OMP) algorithm and and the 1 1 —norm minimization 
(also called basis pursuit) are two widely studied CS reconstruction algorithms [28]. 

Candes and Tao [12] have introduced the following isometry condition on matrices <f> 
and have established its important role in CS. An m x M matrix <f> is said to satisfy the 
Restricted Isometry Property(RIP) of order k with constant 5k (0 < 5k < 1) if for all vectors 
x G M M with ||a;||o < A, we have 

(1 - 5 k ) \\x\\l < || < Pr ||2 < (1 + 5 k ) \\x\\l . (3) 

The following theorem [10] establishes the equivalence between Po and Pi problems: 

Theorem 2.1. Suppose an m x M matrix <f> has the (2k, 5) restricted isometry property 
for some 5 < \[2 — l, then Po and Pi have same A —sparse solution if Po has a A —sparse 
solution. 

The mutual-coherence //(<&) of a given matrix $ is the largest absolute inner-product 
between different normalized columns of <f>, that is, /u( < f>) = max!< jj< m, i^j prfypr^- Here, 
<3>fc stands for the A-th column in <f>. The following proposition [7] relates the RIP constant 
5k and /i. 

Proposition 2.2. Suppose that <3>i,..., &m axe the unit norm columns of the matrix with 
coherence /_/. Then satisfies RIP of order A with constant 5k = (A — T)y. 
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One of the important problems in CS theory deals with constructing CS matrices that 
satisfy the RIP for the largest possible range of A. It is known that the widest possible range 
of k is of the order lo ^My for some well-behaved order constant [5, 15, 24], However the 

only known matrices that satisfy the RIP for this range are based on random constructions 
[5, 10]. Presently the researchers working in CS theory attempt issues [1, 2, 18] such as 

1. Maximizing the sparsity of the solution, that is, for a given pair (y, <3>), hireling x such 
that ||x||o is as small possible 

2. Improving the recovery process, that is, incorporating the prior information (for ex¬ 
ample, partial support solution) into the recovery process 

3. Developing faster algorithms 

4. Designing matrices that satisfy recovery properties 

The recovery algorithms available to date use matrices <f> that are generated at random 
and there is no efficient method for verifying whether or not a given matrix <f> does indeed 
have the stronger reconstruction properties. For this and other reasons, it is useful to have 
an explicit construction of <f\ 

3. Construction Procedure: 

In this section, we present our deterministic construction procedure taking two existing 
binary sensing matrices as inputs. We show that the resulting matrix has smaller density 
than the original two. Let 4 / mX M be a binary sensing block matrix having k" row blocks (the 
meaning of row blocks is explained later with example), each of size n such that every block 
contains single 1 and the overlap between any two columns is atmost r, that is |{/|\Eh jP = 
1} = 1}| — r f° r an y two columns and 4/ g . Let v b( n , xM , be another binary sensing 

block matrix having k! row blocks, each of size n' such that each block contains single 1 and 
the overlap between any two columns is atmost r'. Assume that n' < n , k < min (A:', k"} 
and r < r' < k < n. Now a new matrix can be constructed by the following steps: 

Step-1: Let the i th column of T be T,. For 1 < i < M, define S) = (((supp('F i ) — 1)( 
mod n)) T + 1 T ), where l r is the vector consisting of all ones of size k" x 1, supp( v Ig) is the 
support vector of x I J t . Since T has k" blocks and each block contains one 1, it follows that 
|Si| = k" where S) is a k" — tuple on the set X = {1,2 ,,n}. Since T has M columns, we 
have M such A:"—tuples. For example, suppose m — 9 and T has 3 blocks, then each block 
is of size 3. Now, if the i th column vp* is (1 0 0 0 1 0 0 0 1) T , supp('Ik) is [15 9], then the 
triplet, Si, corresponding to this column is [12 3]. Similarly from the matrix 'L / , we can 
generate M' number of A/—tuples on the set Y = {1,2,..., n'}. 

Step-2 : From A;"—tuples of hrst matrix \P we remove last (k" — k) entries from each tuple 
to make it a A;—tuple. Now add (—1) to each of the entries of the M number of A;—tuples 
that are obtained from T. Let the A;—tuples be of the form (qi, ..., c^) for 1 < i < M. 
From A;'—tuples of second matrix T', we remove last (A/ — k) entries from each tuple to make 
it A:—tuple. By this process, we get M' number of A’—tuples from the second matrix \Effi Let 
each A—tuple have the form (c' 1; c' 2 ,..., c' jk ) for 1 < j < M’. 
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Step-3 : Now replace each k —tuple (cn, Ci 2 , • ■ ■, Cik) with M' number of k —tuples (ob¬ 
tained from T') by adding n'(cn, Ci 2 ,..., Ci k ) to each of the k —tuples (c'- 1; c' 2) ..., c' ]k ) for 
1 < j < M'. This way, we get MM' number of k— tuples on the set X' = {1, 2,..., nn'}. 
Denote these k —tuples by F. Let (ap, aj 2 , ..., a,-*,) for j = 1, 2, MM' be the k— tuples 
in F. 

Step-4 : From these k— tuples we form a binary vector of length knn' where 1 occurs 
in the positions (l — 1 ).nn! + a,ji for l — 1,2,..., k and rest of the positions are zeros. Using 
these MM' number of A;—tuples, we form binary sensing matrix $ having k number of blocks 
and each block is of size nn' and every block contains single 1. The position of one’s in each 
block is indexed by the k— tuples. So the size of the matrix $ becomes nn'k x MM'. The 
pseudo code of the proposed method is as follows: 


1 . 

2 . 

3. 

4. 

5. 

6 . 
7. 


9. 


Input: Two matrices T, T' of sizes nk" x M, n'k' x M' 

Suppose |supp(T,;) fl supp(Tj)| < r and 

|supp(T') fl supp(T')| < r', V i 4 j 

Assume n' > n, k < min{A: / , A:''} and r < r' < k < n 

Set Si = (((supp(Ti) - 1)( mod n)) T + l r ), \Si\ = A;", V i — 1, 2,..., M 

Similarly, S' 3 with |,S'-| = k! are dehned for j = 1, 2,..., M' 

Set S'F = S', u — 1 + n'Si^ki S'j k and Si t k are first k entries of S'' and Si 




1^1 = k and Let ^. = (^,^,...,5; 
1 < S'Fj < nn',\/l = 1 , 2,..., k 


// 


From each S'F, create a vector v % 3 of length nn'k 

and Vij = 1 at ( l — 1 ).nn' + S" 3 l for l — 1 , 2 ,..., k and zero elsewheres 

Output: a matrix of size nn'k x MM', whose columns are v t:J . 


Lemma 3.1. The coherence of the enlarged matrix of size nn'k x MM 1 , yu(4 > ), is atmost 
iL 

k ' 


Proof. Let 4^, 4>j be two arbitrary columns of matrix 4>. There exist two k— tuples /,, f 3 G F 
such that 4>j, are the corresponding vectors of k— tuples /*, fj as defined in the above con¬ 
struction . Suppose fi = (c' fc/1 , 4 , 2 ,..., c' k , k )+n'(c kl , c k2 ,..., c kk ) and f 3 = {d' kn , d' k , 2 ,..., d' k , k )+ 
n'(d k 1 , d k2 ,..., d kk ), where (c' kll , c ' k , 2 ,..., c' k , k ), (d' k , 1 , d ' k , 2 ,..., d' k , k ) are two A:—tuples obtained 
from T' and {c k \, c k2 ,..,, c kk ), (d k i, d k2 ,..,, d kk ) are A:—tuples obtained from T. Set g tJ = 

fi~ fj = ( c fe'l5 c 'k'2i ■ ■ ■ 1 C 'k'k) — d'k’2i ■ ■ • ) ^k'k) + n ' ^( C fcl 5 C k2 , • • • > C kk ) — (d k i, d k2 , . . . , d kk )^j . 

We will show that g l3 has atmost r' number of zero elements, which implies that intersection 
between the supports of 4>j,4>j is atmost r'. Now \c',. — d',.\ < n' and \c ki — d ki \ < n. The 
I th element of g t j is (gij)i = ( c' kn — d' k u ) + n'(c k i — d k i ). We investigate the cases wherein 
G 9ij)i = 0- 

Case 1 : Suppose C 'k'l 4 d 'k'i- 

Sub case 1.1: If c ki = d kh then (g l3 )i 4 0. 
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Sub case 1.2: If cj~i fi du, then also (gij)i 0. Since 1 < \cki~dki\ < n and \c ',— d ',| < n', 
we have (g^i = (c^ - + n'(c w - 4z) 7^ 0. 

Case 2: Suppose c' ¥l = d! kn then (gij)i = 0 only when cu = du ■ 

From the above cases, we conclude that (gij)i = 0 only when c' ¥l = d ¥l and Cki = dki, 
but as both matrices T and T' can have atmost r' intersections, (gij)i = 0 can occur for 
atrnost r'— choices of j. So have atmost r'— intersections. So the coherence yn( < h) of <h 

is atmost j. 

□ 

The following theorem shows the RIP compliance of 

Theorem 3.2. The afore-constructed matrix nn'kxMM> satisfies RIP with 4' = (k 1 — 1)( j) 
for any k' < p + 1. 

Proof. Proof follows from the Proposition 2.2 and Lemma 3.1. □ 

Remark-1: The density of the matrix is -4 which is smaller than the densities of the 
matrices T and T', which are ^ and ^ respectively. Hence using our construction procedure 
we can construct more sparse binary matrices from existing binary sensing matrices. 
Remark-2: The positions of l's in the matrix are known and we know that every block 
contains single one. So we can permute the columns of the matrix so that the support of 
the permuted matrix has nice structured sparsity. 

3.1. Example: 

In this section, we demonstrate our proposed method of construction via an example. 
Let ^4x4 be a binary sensing block matrix 


/ 1 0 1 0 \ 

0 10 1 

10 0 1 

\ 0 1 1 0 / 

having 2 row blocks, each of atmost size 2 such that every block contains one 1 and the 
overlap between any two columns is 1. Let \Hg xg be a binary sensing block matrix 


/ 100100100 X 
0 1 0 0 1 0 0 1 0 

001001001 
1 0 0 0 0 1 0 1 0 

0 1 0 1 0 0 0 0 1 

001010100 
1 0 0 0 1 0 0 0 1 

010001100 
V 001100010 / 



having 3 row blocks, each of atmost size 3 such that every block contains one 1 and the 
overlap between any two columns is 1. 

Step-1 : The corresponding four 2—tuples obtained from T constitute S as 


5 = 



The corresponding nine 3—tuples obtained from T' constitute S' as 



Step-2 : Now removing the last entries from each 3—tuple of S', we get S" that contains 
nine 2—tuples as 


S" 







3 

2 


Step-3 : Let us add (—1) to each of the entries of S and then multiply each of the tuples 
by 3. Now to each 2—tuple add all 2—tuples of S" and generate S"' with IS""! = 36. Now 
all the entrices of S'" are less than or equal to 6. All the tuples of S'" are as follows: 


J fl\ /4\ fl\ (4\ (2\ (5\ (2\ (5\ (3\ ( 6\ f3\ ( 6\ /1\ /4\ /l\ /4\ f2\ f5\ (2\ ( 5\ f3\ (6\ (3\ ( 6\ 

V 4 / V 4 /’ V 1 /’ V 2 / V 5 /’ V 5 /’ V 2 /’ V 3 /’ V 6 / V 6 /’ V 3 /’ w’ V 5 /’ V 5 /’ \ 2 /’ V 3 / V 6 /’ V 6 / V 3 /’ V 1 /’ V 4 / V 4 / V 1 /’ 

fl\ f 4\ fl\ (4\ (2\ (5\ (2\ (5\ (3\ f 6\ (3\ f 6\ 1 

1^3 J ’ y6 J ’ y6 J ’ \^3 J ’ ^4 J ’ \^4 J ’ \^2 J ’ ^5^ ’ ^5^’ ^2^ j ’ 


Step-4 : From each 2—tuple of S'”, we form a column consisting of two blocks of size 6 
where the positions of l's in each block are dictated by the two entries of the tuple. Finally, 
from all these steps, we obtain the following matrix $12x36 : 


/10100000000010100 
00001010000000001 
00000000101000000 
01010000000010100 
00000101000000000 
00000000010100000 
10010000000000000 
00001001000010010 
00000000100100001 
01100000000000000 
00000110000001100 
\0000000001 1000000 


00000001010000000 

01000000000101000 

00010100000000010 

00000000101000000 

10100000000010100 

00001010000000001 

00010010000100100 

00000000000000010 

00100001001000000 

00001100000011000 

00000000000000001 

11000000110000000 


0 0 \ 
0 0 
1 0 
0 0 
0 0 
0 1 
0 0 
0 1 
0 0 
0 0 
1 0 
0 0/ 


3.2. Example 

Let pi,2 be two distinct primes and choose a positive integer r, such that r < p 2 < Pi- 
Using the method in [15], one obtains binary matrices T, T' of sizes p\ x p r 2 +1 and p\ x p\ +1 
respectively. Thus we get p\ +1 number of pi—tuples and p 2 +1 number of p 2 —tuples. If 
we apply our construction procedure on these matrices then we generate matrix d> of size 
Pip| x (pip 2 ) r+1 with coherence The density of this matrix $ is which is small 
compared to — and —, the densities of T and T' respectively. 

In the construction of 4>, the row size is P1P2 and each column contains p 2 number of ones, 
and the overlap between any two columns is r , therefore the maximum possible column size 

(P1P2\ (PIP2\ 

for this construction [4] is , which is of the order (PiP2) r+1 , that is, = 0((pip 2 ) r+1 ). 

L+lJ Vr+lJ 

Consequently, our enlarged matrix attains the maximum possible column size asymptotically. 
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4. An application of proposed construction procedure: 

As stated already, the present work helps us to construct deterministic CS matrices of 
lower density. In addition, we show that this methodology could be adopted to constructing 
CS matrices of general row sizes. 

Theorem 4.1. Suppose m is any positive integer different from p,p 2 ,pq for distinct primes 
p, q. Then there exists a binary CS matrix of row size m. 

Proof. Case 1: If rn — //, i > 2, then m can be written as m = p l ~ 2 p ■ p (let take p = k < 
n' = p < p l ~ 2 = n). Suppose ft is a matrix having p'(> p)— blocks where each block is of size 
p* -2 with intersection between any two different columns is r. Let M be a CS matrix having 
p blocks where each block is of size p with intersection between any two different columns is 
r'(It is guaranteed in [15] that such matrices exist and in above Example 3.2). By applying 
our construction procedure on these two matrices, we generate matrix <f> of row size p*. 
Case 2: Suppose m = pfpf such that both r\ 1 r 2 and p[ l > pf. Take a CS matrix T 
having pi— blocks with each block being of size pf~ l - Let M be another CS matrix having 
Pi—blocks where each block is of size pf (existence of such matrices is guaranteed in [15]). 
By applying our construction procedure on these two matrices, we generate matrix <f> of row 
size rri. 

Case 3: Suppose rn = p'f pf • • • p(“ such that s > 3 and pf > pf > ■ ■ ■ > p r s a . By induction 
on s, we show that there exists a CS matrix of row size m. For s = 3, set m 3 = pfpfpf. 
Suppose the CS matrix T has pg 3 —blocks where each block is of size p[ l . Let M be another 
CS matrix having pi] 3 —blocks such that each block is of size p r f (existence of such matrices 
is guaranteed in [15]). If we apply our construction procedure on these two matrices, then 
we generate matrix <f> of row size pfp'fpf , which is m 3 . Assume that for 3 < l < k, there 
exists a matrix of row size mi = pfpf ■ ■ ■ pf ■ Now we have to show that there exists a 
CS matrix of row size rriy. = pf pf ■ ■ ■ pflfVkClP^ ■ Take a CS matrix T having p r f —blocks 
such that each block is of size m.k -2 = pfpf '' ' Vk^l ■ Suppose T' has pff —blocks such that 
each block is of size pfl{. By applying our construction procedure on these two matrices, 
we generate matrix <f> whose row size is pf pf ■ ■ ■ p^flfVkClVk 1 which is m. By induction we 
conclude that there exists a CS matrix of row size m = pfpf ■ ■ ■ p r s a . □ 

Remark-3: In the process of obtaining a binary matrix of row size m if we use matrices 
Ti,..., Ti with column sizes Mi,..., M t respectively and with coherence pi,..., p f respec¬ 
tively, then the resulting matrix $ is of size m x (Mi • • • M t ) and the coherence p$ of $ is 
max{pi,... ,p t }. 

Remark-4: The authors of [3] obtained general row size via Kronecker Delta product. If we 
apply Kronecker Delta product on \l I mX M and T, the existing matrices that we started 
with, one obtains T", which is of the size mm' x MM'. The method proposed herein, on 
the other-hand, results in <f> of size nri'k x MM'. Although both the matrices have same 
coherence and same density, as nri'k < mm', has better aspect ratio or redundancy factor 
(column to row ratio) compared to M. 

Remark-5: As coherence of the resulting matrix depends on the choice of the two matrices, 
we can choose the two matrices in such a way that the resultant matrix has low coherence. In 



the above proof, we have used as an example the matrices constructed in [15]. Nevertheless, 
one may use other suitable binary matrices like those constructed in [23]. 

5. Construction of ternary matrices 

Let $ mX M be a matrix having k number of blocks where each block is of size n, containing 
single one’s and intersection between any two different columns of $ is atrnost r. Let the i th 
column of <f> be For 1 < i < M, define f = (((supp(<f>j) — 1)( mod n)) T + 1 T ), where 
1 T is the vector consisting of all ones of size k x 1, supp( < hj) is the support vector of 
Since <f> has k blocks and each block contains one 1, it follows that If = k where f t is a 
k— tuple on the set X = {1,2,..., n}. Since <L has M columns, we have M such k— tuples. 
Let fj = (f n . fi 21 ..., f ik ) be the k— tuples. Now decompose this k— tuple into k —number 
of two tuples as (l,fi t ), where l denotes the I th block position and f denotes the I th entry 
of fi. It is to be noted here that 1 < l < k, 1 < f < n. 

5.1. Construction 

From the k— tuple f = ( f tl . f i21 , f lk ), we form a binary vector of length n x k where 
1 occurrs in the positions (l — 1 ).n + f for Z = 1,2,...., k and rest of the positions are zeros 
(the vector is nothing but ) and then we replace 1 with (—1) if l > f H , which results in 
a vector consisting of 0,1, —1 as entries. As we have M such k— tuples, we get M number 
of ternary vectors of length n x k. This way, we get a ternary matrix < h( nxM of size same as 
that of $ mxM - 

Let (0 1000110 Of be a column from a matrix having 3 blocks and each block is of 
size 3. Then the support of this column is the 3— tuple (2 3 1). Break this 3- tuple into three 
2- tuples (1 2), (2 3), (3 1) as we mentioned above. Now in the 3rd block as 1 occurs in the 
1st position, we replace 1 with —1 in the 3rd block of the binary vector (0 1000110 Of 
to generate the ternary vector (0 1 0 0 0 1 (—1) 0 Of. 

Remark-6: The coherence and density of the ternary matrix remain same as that of the 
binary matrix <L. 

As in the case of binary CS matrices, it is possible to construct ternary matrices of 
general row size as well, which is concluded as the following theorem: 

Theorem 5.1. Suppose m is any positive integer different from p, p 2 , pq for distinct primes 
p, q. Then there exists a ternary CS matrix of row size m. 

Proof. Initially we construct a binary matrix of row size m, which is discussed in Section 4. 
Then we use methodology in section 5.1 to obtain a ternary matrix of row size m. □ 

5.2. A different approach to constructing ternary matrix: 

In this section, we present another deterministic construction procedure of ternary matrix 
by combining binary and Hadamard matrices . We show that the resulting matrix has same 
density and coherence as binary but with better aspect ratio (column to row ratio). Let 
v I / mx m be a binary sensing matrix having k number of one’s in each column and overlap 

9 



between any two columns is r. Suppose there exists a Hadamard matrix H of size (k + r') 
for some r' e {0,..., r}. For each column of T, we replace each of its 1-valued entries with 
a distinct row of H to obtain a ternary CS matrix <F of size m x M{k + r'). As the rows of 
Hadamard matrix H are orthogonal, the rows of $ are orthogonal. From the construction 
methodology, it is easy to check that the coherence of the matrix $ is £. The density of the 
matrix $ is A. 

m 

6. Concluding Remarks: 

As CS matrices of low density (or high sparsity) result in algorithms with low compu¬ 
tational complexity, the present work has constructed sparse CS matrices from the existing 
two binary sensing matrices. As an application of proposed methodology, we have shown 
that binary and ternary CS matrices for a more general set of numbers can be constructed. 
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